PROCEEDINGS 

°F SCIENCE 



Fade approximants and g - 2 for the muon 



T— I 
O 
(N 
> 

o 

O: 
a^: 

(N 



Christopher Aubin 

Department of Physics and Engineering Physics, Fordham University 
Bronx, NY 10458, USA 

Thomas Blum 

Physics Department, University of Connecticut 
Storrs, CT 06269, USA 

Maarten Golterman*^ 

Institut de Fisica d'Altes Energies, Universitat Autdnoma de Barcelona, 
E-08193 Bellaterra, Barcelona, Spain 



l> 
d 

(N 



Santiago Peris 

Department of Physics, Universitat Autdnoma de Barcelona 
E-08193 Bellaterra, Barcelona, Spain 



The leading hadronic contribution to the muon anomalous magnetic moment is given by a 
weighted euclidean momentum integral of the hadronic vacuum polarization. This integral is 
dominated by momenta of order the muon mass. Since in lattice QCD it is difficult to compute 
the vacuum polarization at a large number of low momenta, a parametrization of the vacuum 
polarization is required to extrapolate the data. Most fits to date are based on vector meson domi- 
nance, which introduces model dependence into the lattice computation of the magnetic moment. 
Here we introduce a model-independent extrapolation method, and present a few first tests of this 
new method. 
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Figure 1: Left panel: the weight f{Q^) for the muon; right panel: typical data for Tl{Q^) from the lattice. 
Horizontal axis: in GeV^. 



1. Introduction 

The anomalous magnetic moment of the muon g — 2 has been measured with great accuracy 
[|l|], and will be measured with even greater accuracy in the near future. Therefore, a reliable 
computation of = g — 2 from theory with a comparable error would provide a precision test of 
the Standard Model that is sensitive to a large class of models of new physics beyond the Standard 
Model. For this reason there has recently been a lot of interest in lattice computations of with 
controlled errors ^, ^, ^; for an overview, and more references, we refer to Ref. [^. Here, we 
report on recent work on the leading hadronic contribution to , which comes from the hadronic 
vacuum polarization [^. 

The contribution to from the lowest-order hadronic vacuum polarization can be written as 
an integral over the subtracted vacuum polarization n{Q^) - n(0) as a function of euclidean Q 

af^'' = Aa'j^ JGV(e')(n(o)-n(e')) , (1.1) 

where f{Q^) is a kinematic weight shown in Fig. |j, left panel. The right panel shows a typical 
example of lattice data for Yl{Q^) (these are the data from a x 144 lattice with lattice spacing 
0.06 fm and niji = 220 MeV discussed below). 

Clearly, one needs to fit these data in order to compute the integral. In most lattice computa- 
tions of a^^^ to date this has been done with various variants of vector-meson dominance (VMD).^ 
This introduces model dependence into the computation, and the aim of the work presented here is 
to remove this model dependence. 

2. Multi-point Pade approximants 

We start from the observation that we can write a subtracted dispersion relation for n(Q^): 

(n(o)-n(e^))/e==£,.,^.*(e^), (2.1) 

' Ref. 1^ used Pade approximants, but, as we will see below, of a different type than those supported by a conver- 
gence theorem. 
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in which p (t) = lmll{t) / K is the spectral function. Because the spectral function p{t) > 0, the 
integral in Eq. is a Stieltjes function, analytic everywhere except along the cut {—°°, — 4m^]. 



For such a function, there exists a theorem, proven in Refs. [ |10| , |11|]: 



Theorem: Given P points (2^,<I>(G- )), / G {1, . . . a sequence of Pade approximants (PAs) 
can be constructed which converge to ^{Q^) on any closed, bounded region of the complex plane 
excluding the cut, in the limit P — oo. 



This sequence of PAs can be constructed from the P points through a continued fraction 

(!2^-6f)^i(6l) 



^iQ') = .J^SL .... , (2.2) 



1 + 



1+ _ 

■ (e^-6^_i)'t'f-i(e^) 



with ^'i related to <I>(ej<;+i) {^q = ^(Qf), etc.). Equation yields a [[{P - \)/2],[P/2]] 
PA (where [x] is the integer part of x). Furthermore, one can prove that this can be rewritten as 

n(e2) = n(o)-G2U+I , (2.3) 

with 

a„ > , bip/2] > ...>bi>4ml , (2.4) 

i.e., all poles are single poles, they are located on the cut, and all residues are positive. The constant 
ao = for P even. 

In the situation of an actual fit to data for Tl{Q^) obtained from a numerical computation, these 
data are only known within some statistical errors. That implies that we do not know any points 
of the function exactly, and a multi-point sequence of PAs as implied by the theorem cannot be 
constructed. Our strategy will be to fit a fixed number of data points on a given interval, using 
the fact that since Il{Q^), according to the theorem, can be described by a converging sequence of 
PAs of the form (pj|), this equation provides a valid functional form to which to fit the data. More 



concretely, we will fit the form (^) for P G {2,3,4,5}; this yields [0, 1], [1, 1], [1,2] and [2,2] PAs. 
In order to compare diffferent fits, we will then compute 

-1 GeV- 



4a^j^ dQ^fiQ^) {n{0) -n{Q^)) . (2.5) 

We note that VMD is the same as a [0, 1] PA, but keeping bi = fixed: This is not a valid PA in 
the sense of the theorem, because the theorem does not say anything about the possible values of 
the parameters in addition to the conditions (^!4|). 

3. Tests 

For our first test, we explore fits to a MILC data set on a 28^ x 96 lattice with lattice spacing 
0.09 fm, and a pion mass m| ^ 480 MeV 1 13]. This is one of the data sets that was also used in 
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correlated 


uncorrelated 






interval < < 0.6 GeV^ 


interval < 2^ < i GeV^ 


PA 


# parameters 








10i«a«L°'2^^i 


VMD 


2 


5.86/3* 


363(7) 


4.37/18 


413(8) 


[0,1] 


3 


11.4/8 


338(6) 


3.58/17 


373(37) 


[1,1] 


4 


7.49/7 


350(8) 


3.36/16 


424(116) 


[1,2] 


5 


7.49/6 


350(8) 


3.35/15 


443(293) 


[2,2] 


6 


7.49/5 


350(7) 


3.35/14 


445(432) 



Table 1: Results for various fits. The fit marked * was done on an interval < < 0.35 GeV^. 




q' (GeV' ) q' (GeV' ) 



Figure 2: Correlated (solid curve) and uncorrelated (dashed curve) for the [1,1] PA fits (left panel), and for 
the correlated [1,1] (solid curve) and uncorrlated VMD (dashed curve) fits (right panel). 

Ref. [|]. Weshow the results for thePG {2,3,4,5} PAs and for VMD in Table The uncorrelated 
VMD fit is the same as the fit to these data performed in Ref. and the results agree. 
Table [I] leads us to make the following observations: 

• The conelated VMD fit is a bad fit as measured by per degree of freedom (dof); adding 
parameters the fits clearly improve. Note that we always choose the fitting interval by looking 
for a minimal value of ;^^/dof. 

• It turns out that it is difficult to determine the parameters of the second pole with any precision 
[^] (as can be inferred from the values of ;^^/dof), but a^^'^ is insensitive to the second 
and higher poles. 

• There is good internal consistency between all fits shown in the table, except between the 
uncorrelated VMD fit and any of the correlated PA fits. However, the VMD fits are model 
dependent, which translates into an unknown systematic error in these fits. 

We display some of the fits of Table |I] in Fig. ^. Not surprisingly, the uncorrelated fits look 
better at small Q^, but all fits shown in the figure do a good job of describing the data. Therefore, 
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Figure 3: The integrand of Eq. (1.1), using the correlated [1,1] PA fit to the 64^ x 144 data set (solid curve), 
compared with the data weig hted by/(e2) in Eq. (1.1). 



based on the data, it is not possible to decide which of these fits is the best fit. 

We repeated our explorative analysis on MILC lattices on a 64^ x 144 lattice with lattice spac- 
ing 0.06 fm, and nijc ~ 220 MeV. We find very similar results;^ in particular we find 

^hlo,q2<i ^ 572(41) X 10"'° , [1,1] correlated , (3.1) 
^hlo,q2<i ^ g4gj.g^ ^ ^Q-io ^ yjyjQ uncorrelated . 

Our conclusions are the same as before. We note that for both data sets the discrepancy between 
the correlated [1, 1] PA and the uncorrelated VMD fit is about 15%. From the point of view that 
both types of fit give a good description of the data, we take this to imply that there is a systematic 
error of (at least) this size afflicting the determination of a^^^ from the lattice. 

The underlying problem is displayed in Fig. ||, where we see that there are essentially no data 
in the region dominating the integral in Eq. ( |1.1[ ). For this, one clearly needs data at more low 
values of Q^, with smaller errors. It would be interesting to see whether these improvements can 
be attained by using twisted boundary conditions, something that has been tried in this context in 
Ref. [U, and by an error reduction technique such as that proposed in Ref. []14|]. 



We have also compared our PA fits with polynomial fits; results are shown in Table ^. "Poly n" 
indicates a fit with a polynomial of degree n. All fits are correlated fits; and the pairs of fits "Poly 3," 
"[1, 1]" respectively "Poly 4," "[1,2]" have the same number of parameters. We observe that the fits 
deteriorate in the polynomial case going from Poly 3 to Poly 4, with errors increasing, and central 
values for a^^^"^ fluctuating more, while this is not the case going from the [1, 1] to the [1,2] 
PA fit. 



^Of course, central values of a^^'~''^ are quite different, if only because of the smaller pion mass. 
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Poly 3 


Poly 4 


PA [1,1] 


PA [1,2] 


# points 










X^/dof 


4^ 


X^/dof 


4^ 


16 


9.6/12 


543(35) 


9.5/11 


483(244) 


9.7/12 


564(55) 


9 Jin 


565(41) 


18 


11.4/14 


526(33) 


10.5/13 


596(79) 


11.2/14 


541(46) 


11.5/13 


561(21) 


20 


13.1/16 


536(23) 


13.1/15 


535(45) 


13.9/16 


572(41) 


13.9/15 


572(37) 


22 


16.5/18 


541(23) 


15.9/17 


513(44) 


18.5/18 


566(37) 


18.5/17 


566(33) 


24 


16.6/20 


537(18) 


16.4/19 


521(41) 


19.4/20 


583(34) 


19.4/19 


583(33) 


26 


30.7/22 


505(16) 


23.6/21 


580(32) 


26.8/22 


557(31) 


26.7/21 


560(27) 



Table 2: Comparison of polynomial with PA fits, abbreviating ^ by . The number of data 

points included in the fit is indicated in the first column, with 20 points corresponding to the fitting interval 
< < 0.53 GeV^. Data from the 64^ x 144 MILC lattices. 



4. Conclusions 

We presented a new method for parametrizing the momentum dependence of the hadronic 
vacuum polarization, with the aim to avoid the model dependence of the VMD-based fits that up to 
this point have been used in most fits to lattice data for the hadronic vacuum polarization. It turns 
out that this is possible, because the vacuum polarization can be represented in terms of a Stieltjes 
function, for which sequences of Pade approximants can be constructed which converge uniformly 
to the function on any bounded region in the complex place excluding the cut. 

We have tested this new idea on two examples of lattice data for the vacuum polarization. We 
note that the fits based on Pade approximants can lead to larger statistical errors than some of the 
VMD fits, as for instance in Eq. (^). However, it should be emphasized that the latter are afflicted 
with an unknown systematic error originating in the inherent model dependence of VMD-based 
fits. The fits based on Pade approximants avoid this systematic error.^ 

The new method looks promising. However, it is clear that data for the hadronic vacuum 
polarization at more low values (of order the square of the muon mass), and with smaller errors, 
will be needed in order to reach a higher precision for a^'^ . As we have seen, fits based on Pade 
approximants and VMD-based fits (both correlated and uncorrelated) give a good description of 

HLO *C 1 

the data, but lead to values for - which differ by about 15%. 

Finally, we observe that ^ — 2 is an example of a quantity which is quite sensitive to the value 
of the pion mass. Therefore, better data for the hadronic vacuum polarization will also have to be 
obtained at small values of the pion mass, certainly significantly smaller than 300 MeV 
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